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Theory of magnetotunneling spectroscopy in spin triplet p-wave superconductors 
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We study the influence of a magnetic field H on the zero-bias conductance peak (ZBCP) 
due to zero-energy Andreev bound state (ZES) in normal metal / unconventional supercon- 
ductor. For p-wave junctions, ZBCP does not split into two by H even for sufficiently low 
transparent junctions, where ZBCP clearly splits for d-wave. This unique property originates 
from the fact that for p-wave superconductors, perpendicularly injected quasiparticle form ZES, 
which contribute most dominantly on the tunneling conductance. In addition, we show that 
for pa;-|-ipy-wave superconductor junctions, the height of ZBCP is sensitive to H due to the 
formation of broken time reversal symmetry state. We propose that tunneling spectroscopy in 
the presence of magnetic field, i.e., magnetotunneling, is an promising method to determine 
the pairing symmetry of unconventional superconductors. 
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Nowadays, the formation of zero-energy Andreev 
bound states (ZES) at surfaces""^' or interfaces^-* of un- 
conventional superconductors is receiving increasing at- 
tention. This state, which is created by injected and 
reflected quasiparticle's feeling different signs of the pair 
potential, can play an important role in determining the 
pairing symmetry of anisotropic superconductors. The 
detection of the ZES is reflected as an observation of a 
zero-bias conductance peak (ZBCP) in tunneling conduc- 
tance, which is generally considered as a clear signature 
of anisotropic pairing.*' ^"^^ The ZBCP is observed in 
various superconductors having anisotropic pairing sym- 
metry, such as the cuprates.^' Sr2Ru04,^'*' and 
UBeia.^^-* Possibility of ZBCP has also been theoreti- 
cally predicted for organic superconductors (TMTSF)2Ar 
very recently. ^^"^^^ The existence of the ABS has crucial 
influences on many transport phenomena in unconven- 
tional superconductor junctions.'^' 

Since the existence of ZES is a universal phenomena 
expected for unconventional superconductors having pair 
potential that changes sign on the Fermi surface, dif- 
ficulty may arise in determining the pairing symmetry 
only from conventional tunneling spectroscopy. 

In order to overcome this difficulty, we require some 
in situ way of probing the symmetry from the tunneling 
spectroscopy. Here, we will show that a promising way 
is to use a tunneling spectroscopy in the presence of a 
magnetic field. For (ia;2_j^2-wave junction, it has been 
shown that screening currents shift the ZES spectrum 
(Doppler shift) and lead to a sphtting of ZBCP. 26-29) 

By contrast, we show in the present paper that for p- 
wave cases, ZBCP does not split into two in the presence 
of a magnetic field since the most dominant contribution 
to tunneling conductance is given by perpendicular injec- 
tion, where the energy shift of the quasiparticles does not 
occur at all because the component of the Fermi velocity 
parallel to the interface is zero for a cylindrical Fermi 
surface. We also show phenomenologically the absence 



of Doppler shift for more general shapes of the Fermi 
surface, assuming inversion symmetry of crystal, where 
the component of the Fermi velocity parallel to the inter- 
face always have the same magnitude but different signs 
between perpendicularly injected and reflected quasipar- 
ticle states. Finally, we show that for a pj;-|-ipj^-wave 
superconductor, the magnitude of ZBCP is enhanced or 
suppressed depending on which way (-|-z or — z direction) 
to apply the magnetic held. We propose that this depen- 
dence on the magnetic field can be used in detecting the 
broken time reversal symmetry superconducting state. 

We calculate tunneling conductance in normal metal 
/ unconventional superconductor junctions by solving 
Bogoliubov-de Gennes (BdG) equation using quasiclas- 
sical approximation as in our previous theory.'^' As 
regards the triplet pairing cases, we assume that Cooper 
pairs are formed by electrons having antiparallel spins, 
but an extension to more general cases including parallel 
spin pairing or non- unitary cases^'^' is straightforward. 
Now, we consider the case where a specularly reflecting 
surface or interface run along the y-direction. The insu- 
lator located at the interface between normal metal and 
c?-wave superconductor is expressed using U{x) as 



U{x) 




-d, < X <0 
x>0 



where the width and height of the barrier are di and Uq, 
respectively. The magnetic field is applied parallel to 
the z-axis, so that the vector potential can be chosen as 
A{r) = (0, Ay{x), 0). The pair potential depends only on 
X since the system is homogeneous along the y-direction. 
Since we are now considering the situation where the co- 
herence length of the pair potential ^ is much smaller 
than the penetration depth of the magnetic field A, we 
can ignore the spatial dependence of Ay{x) in the fol- 
lowing calculations. We assume Ay{x) = Aq = —HX, 
where H is the applied magnetic field. The normalized 



2 



Letter 



tunneling conductance aT{eV) under the bias voltage V 
is given by^' 



<tt(-E) 
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' (1 - Z2)2 sinli2(Aodi) + 4^2 cosh2(Ao(ii) ' 

, Ao = \/2mUoln^, K, = Kf/Xo, 



\/l — k;2 cos2 



l + aN(0)|r+|2 + [aN(^) 



|i + [c7NW-i]r+r_|2 
MM E±>o, 

-E±<0, 



El-\Aof{e±) |2 + ^± 
-Ao/(0±) 



(3) 



(4) 



^i- I Ao/(0±) |2 



^^4 
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with the injected and reflected angles being 9+ = 9 
and e. = n-9, respectively, Hq ^ ^ = 
^0 = ^i'F/(7rAo), -B = eV, and 5± = vpy/vp, where 
k-p and up is the Fermi momentum and the Fermi veloc- 
ity, respectively, taken to be common in superconductor 
and normal metal for simplicity. Aq is the magnitude of 
the maximiim value of the pair potential and E± are the 
Doppler shifted energies of injected and reflected quasi- 
particle states. In the above, we have assumed a spatially 
constant pair potential in the superconductor. As shown 
in our previous paper, if we concentrate on the qualita- 
tive feature of (Tt(-E') at low voltage, this assumption is 
justified.^' ^"^^ 

Wc take Xodi = 3 and k = 0.5, where the transmission 
coefficient perpendicular to the interface is about 0.02. 
In the following, we will look at aT{eV) and as{9,eV) 
for various cases. 

We have calculated aT{eV) for a da,2_y2-wave junction 
with (110) oriented surface and for a p^^-wave junction 
with (100) oriented surface, where f{9) is expressed as 
f{9) = sin(2^) and f{9) = cos{9), respectively, as a pro- 
totype of d-wave and p-wave superconductor junctions. 
In dx2-y2-wa.ve junctions, aT{eV) has a peak splitting 
[see curve b and c in Fig. 1(a)] in the presence of magnetic 
field iJ29, 30) gjjicg lY^Q transparency of the junction is suf- 
ficiently small. On the other hand, in the p^^-wave junc- 
tions, there is no ZBCP splitting with the increase of the 
magnitude of H, where only the height of the ZBCP is 
reduced and the width becomes broad [Fig. 1(b)]. Thus, 
the response of ZBCP to the magnetic field is quite dif- 
ferent from that in the d.r2_,,2-wave case. 
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Fig. 1. The normalized tunneling conductance aT(eV) in the (a) 
d^2_y2-wavc state with (110) oriented surface and (b) 

Px - wave 

state with (100) oriented surface. Aorf; = 3 and K = 0.5. a: 
H = 0,h: H = 0.2Ho, and c: H = OAHq. 
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Fig. 2. Angle resolved tunneling conductance in the superconduct- 
ing state <7s{0,eV) for = with Xgdi = 3 and k = 0.5. (a) 
d^2_y2-ws,ve and and (b) px-wave. a: = — 7r/4, h:0 = — 7r/8, 
c: 6» = 0, d: 6» = tt/S and e: 6» = 7r/4. 



In order to tmdcrstand this difference, we look into 
C7s(^, eV), i.e., the angle resolved conductance in the su- 
perconducting state. For H = 0, the height of ZBCP 
as{9,0) is given as (Ts{9,0) = 2 cos 6* for f{9) ^ in- 
dependent of a-f^{9) due to the formation of ZES and 
£7s(6',0) = <7iv(6') for f{9) = 0, respectively.^) In the ab- 
sence of the magnetic field, as{9,eV) always has ZBCP 
except for special cases, i.e., 9 = for the dx2 _y2-wauYe 
junction and 9 = ±7r/2 for both junctions. For the 
dj.2_y2-wavc junction, the width of the ZBCP W takes 
its maximum value for some oblique injection angle 9m, 
while for Px-wave, both W and as {9, 0) becomes largest 
at 6* = 0. 

By applying H, the positions of peak shift from zero. 
In the present case, since positive H is chosen, quasipar- 
ticle energy E increases (decreases) for 6* > (6* < 0), 
so that the peak position of as{9,eV) becomes located 
in the negative (positive) voltage region. However, a 
remarkable feature is that the ZBCP for perpendicular 
injection, i.e. crs(0, eV), is not changed at all by H. This 
is because vpy = at ^ = 0, so that there is no Doppler 
shift for perpendicular injection. For the p^j-wavc case, 
perpendicular injection contributes most dominantly to 
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Fig. 3. Angle resolved tunneling conductance in the superconduct- 
ing state as{8,eV) for H = 0.2H() with Xodi = 3 and K = 0.5. 



(a) (1^2 _y2-wavc and (b) px-wavc. a: 
e = 0, d: e = it/8 and e: = 7r/4. 



-7r/4, b-.e = -it/8, c: 



the low bias behavior of the integrated normahzed tun- 
neling conductance aT{eV), and consequently the ZBCP 
is robust against H. In the abovc\ = at = is 
a consequence of assuming a cylindrical Fermi surface. 
In fact, the absence of Doppler shift for perpendicular 
injection can be shown in more generalized cases with 
VFy{kx,ky = 0) 0. Namely, there is generally an in- 
version symmetry of the crystal, e(fe) = e(— fc), so the 
injected and reflected quasiparticle states have different 
sign of VFy {vFy{kx,0) = —VFy{—kx,0)) for perpendicu- 
lar injection. In order to look into this general situation, 
we consider a phenomenological analogue of aR{6,E) as 



p{E) 



i + a|r+p + [a-i]|r+p|r_ 
|i + [a-i]r+r_|2 



(5) 



Here we take a = 0.1 as a typical value for a low trans- 
parent barrier, and calculate p{E) for two cases, i.e., 
i) VFy{k:o,0) = -VFy{-kx,0) i.e., g+ = -g-, and ii) 
VFy{kx,ky) = VFy{-kx,ky) i.e., g+ = g-. As seen from 
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Fig. 4. Local density of state with (a) ■ypyCfex) = —■!'F3/(~'^a:) 
and (b) vpy{kx) = Vfy{—kx)- a: H = 0, h: H = 0.2Ho, c: 
H = 0.4/fo. 

Fig. 4(a), there is no Doppler shift of ZBCP with the 
increase of H for case (i). The absence of the shift is 
quite different from that of the conventional case (ii) 
[Fig. 4(b)], where the magnitude of the shift of the peak 
position is proportional to the magnetic field. The 
absence of Doppler shift originates from the cancella- 
tion of the additional phase in the product of F+F- due 
to the different sign of E+ and E- at E = 0. The 
cancellation of the additional phase shift of the wave 



function is a quite novel feature originating from the 
difference in the sign of VFy between reflected and in- 
jected quasiparticle. Although the present argument for 
VFy{kx,0) = —VFyi—kx,0) 7^ is pheuomeuological, re- 
cently, we have shown for a tight binding model of a 
possibly triplet superconductor (TMTSF)2X that zero- 
energy peak in the surface density of states does not split 
in the presence of a magnetic field, not only for p-wave 
pairing, but also for a triplet /-wave pairing. 

Finally, we now look into the chiral p-wave {px+^Py- 
wave) case, where f{9) is given as f{9) = cxp{iO). This 
is a pairing symmetry possibly realized in Sr2Ru04. As 
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Fig. 5. The normalized tunneling conductance o'T(eV') in the px + 
ipy-wave state. Xodi = 3 and k = 0.5. a; H = 0, h: H = 0.2Hq, 
c: H = 0.4Ho d: H = -0.2Ho and e: H = -0.4i?o. 

seen in Fig. 5, ZBCP again does not split in the pres- 
ence; of a magnetic field, while its height is enhanced 
for positive H [H in the +z direction; curves b and c), 
but is reduced for negative H [H in the —z direction; 
curves d and e). Such an asymmetric H dependence of 
ZBCP around H = does not appear in dx2_y2-wa.ve 
and Pa;- wave cases. In order to understand this asym- 
metric feature in detail, we also plot in Fig. 6 the an- 
gle resolved tunneling conductance as in Figs. 2 and 3. 
As seen in Fig. 6, the peak position of as{9,eV) has a 
strong dependence. Namely, the p^j+ipj^-wavc pairing 
induces a broken time reversal symmetry state (BTRSS), 
so that the peak positions for ^ 7^ is no longer located 
aX E = O"*^^ even in the absence of a magnetic field. Now, 
in the presence of a positive magnetic field, all the peaks 
for ^ 7^ are shifted toward E = 0, [Fig. 6(b)] while 
they are shifted ainay from E = for a negative mag- 
netic field [Fig. 6(c)]. This unique feature can be inter- 
preted as follows. In p^+ipy state, due to the formation 
of BTRSS, there is a surface current''^) which flows par- 
allel to the interface and the spontaneous magnetic field 
even without applied magnetic field. If the direction of 
the applied magnetic field H is the same (opposite) as 
that of spontaneous field, the effective magnetic field is 
enhanced (reduced) for iJ < (H > 0). Then the result- 
ing magnitude of ZBCP in total conductance a-T{eV) is 
reduced (enhanced) for negative (positive) H. In sum- 
mary, we have calculated tunneling conductance in nor- 
mal metal / unconventional superconductor junctions, 
where dx2_y2-wave, p^-wave and Px+ipjz-wave pairings 
have been chosen as a prototype. We focused on the 
influence of the applied magnetic field on ZBCP. For p- 
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Fig. 6. Angle resolved tunneling conductance in the superconduct- 
ing state <Ts{S,eV) in px+ipy-v/a.ve state with Aorfi = 3 and 
K = 0.5. (a) H = 0, (b) H = 0.2Ho and (c) H = -0.2Ho. a: 
B = -n/4, b: 9 = 0, and c: 9 = 7r/4. 



wave cases, ZBCP does not split into two by magnetic 
field since the most dominant contribution to tunneling 
conductance at zero-energy originates from perpendicu- 
lar injection, where Doppler shift does not occur. The 
absence of Doppler shift has been shown in general situ- 
ations, where the component of the Fermi velocity par- 
allel to the interface has the same magnitude but dif- 
ferent signs between injected and reflected quasiparticle 
states. As for pa;-|-ipy-wave superconductor junctions, we 
have shown that ZBCP does not split, while its height is 
sensitive to the direction of the applied magnetic field, 
which is a consequence of a broken time reversal symme- 
try. This behavior should be observed in the tunneling 
experiments of Sr2Ru04, if it is actually a chiral p-wave 
superconductor . 
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